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Abstract. This contribution intends to give a pedagogical introduction to the topic of dark energy 
(the mysterious agent supposed to drive the observed late time acceleration of the Universe) and 
to various observational tests which require only assumptions on the geometry of the Universe. 
Those tests are the supernovae luminosity, the CMB shift, the direct Hubble data, and the baryon 
■ acoustic oscillations test. An historical overview of Cosmology is followed by some generalities on 

FRW spacetimes (the best large-scale description of the Universe), and then the test themselves are 
discussed. A convenient section on statistical inference is included as well. 
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: INTRODUCTION 

d 

Cosmology is a branch of physics which is experiencing a tremendously fast devel- 
opment lately triggered by the arrival of many new observational data of ever more 
exquisite precision. These findings have been crucial for the improvement in our under- 
standing of the Universe, but the vast amount of knowledge on our Cosmos which we 
have nowadays would never have been possible without the concerted effort by exper- 
imentalists and theoreticians. One of the results of this fantastic intellectual pursue is 
the puzzling discovery that there is in the Universe a manifestation of the repulsive side 
of gravity. This is the topic to which this lectures are devoted, and more specifically I 
wish to address some methods which the community believes are useful for building a 
deeper understanding of why, how and when our universe began to accelerate. Put in 
more modest words, these lectures will dissert about how one can take advantage of 
various observational datasets to describe some basic features of the geometry of the 
Universe with the hope they will reveal us something on the nature of the agent causing 
the observed accelerated expansion. 

As this is quite an advanced topic in Physics (Astronomy), this contribution will build 
tougher as we proceed, so I hope the readers will enjoy to start off from a little historic 
stroll in the science of surveying the skies. For wider historical overviews than the one 
presented here, the two main sources I recommend to you, among the so many available, 
are [1] and [2]. 

Historical records tells us that Astronomy was born basically because of the need by 
agrarian societies to predict seasons and other yearly events and by the esoteric need 
to place humanity in the Universe. This discipline developed in many ancient cultures 
(Egyptian, Chinese, Babylonian, Mayan), but it was only Greek people who cared un- 
derstanding their observations and spreading their knowledge unlike in other cultures. 
Their influence was crucial as the modern scientific attitude of relying in empiricism 
(Aristotlean school) and translating physical phenomena into the mathematical language 



(Pithagorean school) built on their way of approaching science. Unfortunately, Aristo- 
tle was far so influential that his erroneous geocentric view of the Universe was not 
questioned aloud until the XVIth century. 

At the beginning of that century Copernicus started to spread quietly his heliocentric 
cosmological model, and even though he did not make much propaganda about his ideas, 
they deeply influenced other figures, such as Galileo, who introduced telescopes into as- 
tronomy and found solid evidence against the geocentric model. Another very important 
contribution to the subject was made by Kepler, who gave accurate characterization of 
the motions of the planets around the Sun. These findings were later on synthetically 
explained by Newton in his theory of gravitation, which built on Galileo's developments 
on dynamics (the law of inertia basically). 

The XVIIIth and XlXth century brought advances in the understanding of the Uni- 
verse as a whole made of many parts which lie not necessarily in the Solar system. Stars 
began to be regarded as far-away objects in motion, and other objects like nebulae were 
discovered, so the idea of the existence of complicated structures outside the Solar sys- 
tem gained solidity. A discovery very much related to the main topic of this lectures was 
the realization that the combination of the apparent brightness of a star and its distance 
get combined to give its intrinsic brightness; basically the amount of energy reaching 
us in the form of light from a distant source decreases is inversely proportional to the 
square of the distance between the source an us (in an expanding the universe definition 
of distance is not the same as in a static one but this broad way of speaking applies all 
the same) This finding let observers realize that stars were objects very much like the 
Sun, or if you prefer they realized the Sun was nothing but yet another star. 

The most important next breakthrough was Einstein's theory of special relativity, 
which generalized Galileo's relativity to introduce light. On the conceptual realm, this 
was a very revolutionary theory at is warps the notions of space and time, and so it set the 
foundations for the best description of gravity so far: Einstein's theory of General Rela- 
tivity. This second theory was the fruition of Einstein endeavors to unify the interactions 
known to him. In this theoretical framework matter/energy modifies the geometry, and 
in turn geometry tells matter/energy how to move/propagate (paraphrasing Wheeler's 
renowned quotation). Among other predictions this theory made a couple which are 
cornerstones of modern astronomy: gravitational redshift (light gets redder as it moves 
away from massive objects), and gravitational lensing (light gets bent as it passes close 
to massive objects). 

The next important advance come from the side of observations. In 1929 Hubble, 
and after having collected data carefully for almost a decade, presented the surprising 
conclusions that galaxies on average move away from us. This effect is encoded in the 
scale-invariant relation known as Hubble's law: v = Hd, where v is the galaxy's velocity 
and d its distance from us. The positiveness of the quantity H as measured by Hubble 
is precisely what told him the Universe is expanding, this being a discovery which gets 
accommodated nicely in Einstein's theory of general relativity. Interestingly, on learning 
about this finding, Einstein discarded his idea of the necessity of some exotic fluid with 
negative pressure to counteract the attractive effect of usual matter (which would make 
the universe contract). It cannot look but funny from today's perspective that Einstein's 
idea has come to life again, at the end of the day the exotic fluid he imagined is one of 
the possible flavors of what the community calls dark energy [3] these days. 



At the risk of not giving everyone the credit they merit, I will just say that recog- 
nition for the concept of a expanding Universe is due to both theoreticians (de Sitter, 
Friedmann, Lemaitre, Gamow,...) and experimentalists (Slipher, Hubble,..)- However, 
the trampolin for this idea to jump into orthodoxy was put by Penzias and Wilson [4] 
who first detected the cosmic microwave background, and by Dicke, Peebles, Roll and 
Wilkinson [5] who were responsible for the not the less important interpretation of those 
observations. The existence of this radiation and its characteristic black body spectrum 
are a prediction of the Big-Bang theory, it is fair to say that if one combines it with other 
sources of evidence, it is almost impossible to refute it. 

The CMB in an invaluable source of cosmological information (visit [7] for a higly 
recommended site on the subject). It has a temperature of of around 2.73 K, with tiny 
temperature differences of about 10~ 5 between different patches of the sky. These 
anisotropics inform us that the photons forming the CMB where subject to an underlying 
gravitational potential which had fluctuations, and this is just and indication of density 
irregularities that seeded the structure one observes today. On the other hand, the image 
of the features of the CMB in different angular scales strongly favors a universe with no 
spatial curvature (flat or Euclidean on constant time slices) that is, it supports the theory 
of inflation. In addition, the CMB has a saying about the fraction of the different fluids 
filling the Universe. 

Fortunately, the mine of cosmic surprises was far from being exhausted, and it stored 
a diamond of many carats in the form of a discovery which has changed greatly the 
mainstream view of the Universe. Up to the late 1990, no repulsive manifestation 
of gravity had been spotted in Nature. In 1998 astronomical measurements provided 
evidence that gravity can not only push, but pull as well, and the repulsive side of gravity 
got unveiled. Very refined observations of the brightness of distant supernovae [8, 9] 
seemed to hint the presence of a negative pressure component in the Universe which 
would make it accelerate, and then a new revolution started. 

One may wonder at this stage what is the relation of supernovae luminosity and cos- 
mic speed up. Objects in the Universe with a well calibrated intrinsic luminosity (super- 
novae, for instance) can be used to determine distances on cosmological scales. Super- 
novae are very bright objects, and so they result particularly attractive for this purpose 
(they can be 10 9 times more luminous than the Sun so there is hope they will be visible 
from up to perhaps 1000 Mpc [10]). As we anticipated in the last paragraph, in 1998 two 
independent teams reported evidence that some distant supernovae were fainted than ex- 
pected. This involved tracing the expansion history of the Universe by combining mea- 
surements of the recession velocity, apparent brightness and distance estimations. The 
most compelling explanation was (and keeps being as far I am concerned) that their light 
had traveled greater distances than assumed. The orthodox view up to then was that the 
expansion pace of the Universe was barely constant, but supernovae seemed to contra- 
dict this. This unexpected and exciting discovery obliged researchers to broaden their 
mind and accept the Universe is undergoing accelerated expansion. 

I should have been able to have convinced the readers by now that these are exciting 
times to be working on Cosmology, as cosmic speed up is such and intriguing phe- 
nomenon with major open questions such as whether dark energy evolves with time, 
how much of it is there, and if is rather not a manifestation of extradimensional physics. 



The answer to these questions requires cannot be dissociated from the response to a 
perhaps more fundamental question: what is the Universe made of? 

The combination of various astronomical observations tell us our universe is basically 
made of thee major components (see for instance [11]). The most abundant one is dark 
energy [3], so this makes it even more interesting to find out whatever we can about 
it. At the other end the by far least abundant component is baryonic matter, and in the 
middle (as abundance is concerned) we have dark matter, in a proportion comparable to 
that of dark energy. There are various sources of astrophysical giving evidence in favor 
of it. Hints of it is existence are provided by the motion of stars, galaxies and clusters, 
but it is known it also played a crucial role in the amplification of the primordial density 
fluctuations which seeded the large scales structure we observe today and dark matter 
imprints can be found in the CMB as well. Dark matter represents quite a challenge 
as its nature remains a mystery; nevertheless if it were baryonic we know from big 
bang nucleosynthesis helium-4 would get converted to deuterium much easily and CMB 
calculations indicate indicate in addition anisotropics would be much larger so the odds 
are most of its is not baryonic. 

Up to here we have made a very broad introduction to our topic with a little bit of 
history and a little bit of physics, but we must not forget Maths are key to Cosmology 
[12]. We only know how to study the Universe using numbers and equations, but 
of course progress in this direction is done with as many reasonable simplifications 
as possible (if they do not compromise rigor, of course). Einstein equations relate 
geometry and matter/energy content of the Universe. Cosmologists are concerned by 
this relation on a large scales picture so as to understand the expansion of the Universe. 
Those equations are non-linear, so studying them is painstaking unless one exploits the 
observational evidence the regularity of the Universe. Non-linearity of those equations 
makes their study a really hard task so simplifications are a must. The two basic ones are 
that galaxies are homogeneously distributed on galaxies larger than 50 Mpc [13], and 
that the Universe is isotropic around us on angular scales larger than about 10 degrees 
[14]. But this is not enough, those two simplifications, which come from observational 
evidence must be completed with the assumption we occupy no special place in the 
Universe. 

This puts on the track of what we could call the parameterized Universe, cosmologists 
work with a greatly a simplified geometric description of the Universe which emerges 
from the latter assumptions. The models for the sources are of reduced complexity too, 
the most common being perfect fluids and fields with known dynamics. This we have 
on the side of theory, but on the side of observations we have to make our own life 
easier too. When doing observations oriented Cosmology, one uses as a variable the 
redshift z of the electromagnetic radiation received, as it encodes information of how 
much the Universe has expanded between emission and reception. Observational tables 
of geometrical quantities can be given, and then one can test different theoretical values 
of the same quantities corresponding to models of interest. The theoretical predictions 
will depend on the sources assumed, and ultimately it will be possible to estimate the 
suitability of a given dark energy model to observations, but which are the available 
probes of dark energy? 

Basically dark energy can be scrutinized observationally from two main perspectives 
[15]: one possility is doing it through its effect of the growth of structures, another 



one is through its impact on geometrical quantities. We will concentrate on geometrical 
constraints/tests in these lectures as in a way they are those which can perhaps be applied 
with less difficulty, although their simplicity does not mean they are the least interesting, 
for instance, the supernovae test is the only test giving a direct indication of the need of 
a repulsive component in the Universe, whereas the baryon acoustic oscillations test is 
thought to have much information in store. Discussion on these two tests will be given 
later on in these lectures. 

Finally, there is one more direction in which these topic is related to Maths apart from 
the geometrical side of it, statistics plays an important role too. Physics is attractive 
because of its ability to "tame" natural phenomena in the sense that laws of physics bring 
order to the apparent chaos of Nature, Astrophysics is even more attractive because it 
evidences how laws of physics apply outside Earth, which is certainly surprising given 
the manifest differences between Earth and every other locations in the Universe of 
which we are aware. However, since these phenomena occur in places far, far away, and 
sometimes they can only be observed indirectly, there is typically an important degree 
of uncertainty. Thus, research is Astrophysics requires understanding not only Physics, 
but also inference. 

Inferential statistics cares for the identification of patterns in the data taking in account 
the randomness and uncertainties in the observations. In contrast, descriptive statistics 
is concerned with giving a summary of the data either numerically or graphically. Both 
will be needed toward two goals: we need to know optimal ways to extract information 
from the astronomical data, but we also need to know rigorous approaches to compare 
theoretical predictions to observations. 

Our approach will be that of Bayesian inference, and we will justify jut below my 
preference, but it must be admitted there is an old vivid controversy on the definition of 
probability between the two main schools: frequentists and Bayesians. Frequentists use 
a definition based on the possibility of repeating the experiment, but this does not apply 
to the Universe, and this sounds like the reason why the number of Bayesian astronomers 
grows every day, but you should not care for this battle right now; just stick to the idea 
that observations related Cosmology needs to resort to Statistics and that Inference is 
vital for constraining dark energy cosmological models. 

By now you should be able to guess what to expect in the next sections. I will 
present you some basics of FRW (Friedmann-Robertson-Walker) cosmologies, then I 
will devote a great deal of this text to details of geometrical tests, and leave for the last 
but one section a convenient primer on statistics. 

GENERAL RELATIVITY AND FRW 

Einstein's influence in cosmology is paramount, Special Relativity is crucial in high en- 
ergy physics, which plays an important role in astronomy and cosmology, as many pro- 
cesses one studies in these areas are very energetic. But General Relativity is even more 
important for Cosmology as it allows to describe the gravitational interaction, which is 
the one governing the dynamics on planetary, galactic and cosmological scales. Gravity 
in astronomy is mostly treated classically as opposed to quantum mechanically, as no 
successful theory of gravity exist. The situation is different with respect to Special Rela- 



tivity, as the connection with the quantum realm is satisfactorily given by quantum field 
theory. Given that Special Relativity is the conceptual precursor of General Relativity, a 
few lines about it are worth before entering an overview of General Relativity. 

Part of the topics of this section are extensively covered in [16, 17, 18, 19]. 

Special Relativity stands on two pillars. The first one is that the laws of physics are 
the same in all inertial frames, reference systems in which bodies are not subject to 
forces remain at rest or in steady linear motion. The second one is that the speed of 
light in vacuum, c, is the same in all inertial frames as confirmed by the Michelson- 
Morley experiment, but actually anticipated intuitively by Einstein. Inferring conclu- 
sions from those two premises requires treating the quantity ct on the same grounds as 
spatial coordinates (say x : y,z), and this interchangeability of space and time make the 
concept of spacetime emerge. When transforming between inertial frames, the quantity 
ds 2 = c 2 dt 2 — dx 2 — dy 2 — dz 2 remains invariant, ds 2 being the norm of the four- vector 
(cdt, dx, dy, dz). 

Special Relativity requires the laws of physics to be written in terms of four- vectors as 
their norms do not change on going from one inertial frame to another. The quantity ds 2 
is called the line-element, and it quantifies the distance between events of the spacetime. 
The line-element is a quadratic form constructed from a matrix g: ds 2 = gj lv dx^ l dx v for 
jii,V=0,l,2,3In Special Relativity gij = diag(l, — 1, 1,1), but in General Relativity it 
need not be diagonal nor constant. 

This brief account on Special Relativity drives us into the theoretical ground of 
General Relativity. In developing this beautiful theory Einstein was influenced by five 
principles: Mach's principle (geometry or motion do not make sense in an empty 
universe); principle of equivalence (the laws of physics look the same to an observer 
in non-rotating free fall in a gravitational field and to an observer in uniform motion in 
the absence of gravity), principle of general covariance (laws of physics must have 
the same for all observers), correspondence principle (from General Relativity one 
must recover on the one hand gravity Special Relativity when gravity is absent, and one 
must also recover Newtonian gravity when gravitational fields are week and motions are 
slow, actually, this principle reflects the very reasonable need that any new scientific 
framework must be consistent with precursor reliable frameworks in their range of 
validity.) 

According to records, walking the tortuous path from Spatial Relativity to General 
Relativity took Einstein 1 1 years. Let us outline the main elements of the construction: 
gravity can be waived locally and regain Special Relativity, locally gravitational effects 
look like any other inertial effect, test particles are assumed to travel on geodesies 
(null ones if they are photons and timelike ones if they are massive), inertial forces 
are accounted for in the geodesic equations by terms which depend on first derivatives 
of the metric (which plays the role of the potentials of the theory), gravitational fields 
make geodesies converge/diverge as described by some terms in the geodesic deviation 
equation which depend on the Riemann tensor (which depends on second derivatives 
of the metric), all forms of energy act as sources for the gravitational field and this is 
encoded in the Einstein equations. 

In General Relativity tensors play a preeminent role. The Riemann tensor or curvature 
tensor R a bcd determines how geodesic deviates, and upon contraction it gives the curva- 



ture or Ricci tensor Rab — g Rdacb- Further contractions allow deriving the curvature or 
Ricci scalar R = g ab R a b, and finally, the connection between matter/energy and curvature 
is encoded in the Einstein equations 

G ab = SKGT ab /c\ (1) 

where G ab = is the Einstein tensor and T ab is the Einstein tensor. 

At this point we have now enough theoretical machinery to study the dynamics of a 
standard universe, but I must insist on the fact that studying the Universe without mak- 
ing radical but reasonable simplifications would be a intractable problem. Fortunately, 
progress can be made because we are lucky enough to have evidence of its regularity. 
As already mentioned in the introduction, we have evidence on the one hand of the ho- 
mogeneity in the distribution of galaxies on scales larger than 50 Mpc [13], and on the 
one hand CMB experiments inform us on the isotropy around us on angular scales larger 
than a 10 degrees [14]. If one then invokes that our place is the Universe is not special at 
all, then isotropy around all its points in inferred. Finally, there is a theorem in geometry 
which tells us that if every observer sees the same picture of the Universe when looking 
at different directions, then the Universe is homogeneous. 

These assumptions boil down into the (Friedmann-)Robertson-Walker metric. Our 
Universe can be viewed as an expanding, isotropic and homogeneous spacetime, and 
that means its line element reads 

ds 2 = c 2 dt 2 -R(t) 2 (dr 2 + S 2 (r)(de 2 + sin 2 (e)d(j) 2 ), (2) 

with R(t) an arbitrary and S(r) a function which can take three distinct forms. If one 
defines R(t) = a(t)Ro, and then makes r — > t/Rq, a more familiar expression is obtained: 

ds 2 = c 2 dt 2 -a(t) 2 (dr z + S 2 (r)(d9 2 + sin 2 (e)d(i) 2 )), (3) 

where a(t) is a dimensionless quantity (customarily chosen to have value 1 at present). 

The geometry of spatial sections in the FRW metric are also worth some further 
discussion. The function S(r) must be such the spatial sections of the RW geometry have 
constant curvature, so the possibilities are S(r) = {sin{r),r,sinh(r)}. Again, a familiar 
form of the line-element is obtain by making S(r) 2 dr 2 — > dr 2 / ( 1 — kr 2 ), and so the three 
cases are respectively k = 1,0,-1 (as shown in Fig. 1). 

Before proceeding, a little remark is convenient. From here on we will use the so 
called natural units %kG = c = 1 under otherwise stated. Having made this comment in 
passing let us turn our attention to the right hand side of the Einstein equations, i.e to the 
matter/energy content. 

We have seen that simplifications in geometry are required to model the Universe. 
In the same spirit reduction of sophistication in the description of matter/energy is also 
required. Simplicity on the one hand, and consistency with observations in the other, 
suggest adopting the perfect fluid picture (no viscosity is assumed): 

Tab = (P + P)u a u b + Pgab, (4) 

with p, p and u a representing the energy density, pressure and velocity field of the fluid. 
Note that in the rest frame of the fluid 



T ab = diag(p, -p,-p-p). 



(5) 



FIGURE 1. From left to right, geometry of spatial sections in a fiat universe (k=0), a positively curved 
or closed universe (k=-l), and negatively curved or open universe (k=-l) 



Let us now formulate Einstein equations for perfect fluids. There are basically two of 
them in these cases. The first one is the Friedmann equation, 

H 2 = a 1 1 a 1 = £ pi/3 - k/a 2 , (6) 

i 

and it acts as a constraint as d is not free, but is rather subject to the amount of energy 
density and curvature. The second one is the Raychaudhuri equation: 

2(d/a - a 1 1 a 1 ) = - £(p,- + pi) + k/a 2 , (7) 

and it is an evolution equation. The combination of the two equations can be used to 
derive a(t) (or t(a) in the least fortunate cases), once p and p have been specified. 

From Einstein equations one can derive other two important equations. The first one 
is the energy conservation equation 

ptot + 3//(p tot + Aot) =0 p t = Y,Pi Pt = J^Pi- ( 8 ) 

i i 

The second one is the acceleration equation, which tells us about the evolution of the 
spatial separation between geodesies 

2d/a = -(p+3p)/3. (9) 

These preliminaries suggest the interplay between the matter/energy content of the 
Universe and its geometry have a crucial influence in its final fate. If we consider a 
model with matter only then p tot = p m = p /a(Y) 3 , and thus 

a 2 = po/3a(t)-k, (10) 

so, any value of a(t) is consistent for k = 0, — 1 for not for k = 1, there is an upper bound 
to a(t). When a(t) = a cr it = Po/3 the model will begin to collapse (because a > —a/2, 
which signals a maximum in the a versus t plot). In contrast, open or flat universes do 
not experience any particular behavior when the critical density is reached. 



(a) (b) 

FIGURE 2. (a) Regions for the LCDM model, (b) A spherical shell in redshift space. 



The addition of a cosmological constant brings a richer set of possibilities (see Fig. 
2(a)). It is convenient to present the cases using fractional densities: 

Cl k = -k/H 2 a 2 , Q A = A/3// 2 , Cl m = p m /3H 2 , Q m + Q A + Q Jt =l. (11) 

Let us consider now perfect fluid cosmological histories. Many p(p) equations of state 
are considered in the literature, but the linear one, p = wp, stars in popularity. Here you 
have a list some physically meaningful cases in the k = case. 

• Electromagnetic radiation, i.e. photons, w = 1/3, p ~ a~ 4 , a(t) ~ t 1 / 2 

• (Incoherent) matter, aka cosmic dust, w = 0, p ~ a~ 3 , a(t) ~ t 2 ^ 3 

• Vacuum energy, aka cosmological constant, w = — 1, p ~ cons, a(f) ~ 

• Quiessence [20], w = cons < - 1/3 ^ - 1 , p ~ fl- 3 ( 1+w ) , a(f) ~ ? 2 /(3(i+w)) 

There are, of course, more complicated equations of state that have received attention 
in connection with late-time acceleration. A list (with some bias, as that of perfect fluids) 
can be this one: 



• Conventional Chaplygin gas [21], p = —A/p, p ~ y'A + B/a 6 

/ \l/(l+a) 

• Generalized Chaplygin gas [22], p = -A/p a , p ~ ( A + B/a 3 ^ +a > j 

• Inhomogeneous equation of state [23], p = —p —Ap a , p l ~ a ~ A(l — a)loga 

Last but not least, another very popular class of sources in Cosmology is that of scalar 
fields. Usually, they can be interpreted as perfect fluids, but the main difference is that 
the equation of state changes over time. Many scalar field models have been proposed 
(for the early and the late universe): 

• quintessence [24], p = <p/2 + V(0) p = <j} 2 /2-V((j)) 

• k-essence [25], p = V((j)) ^2<p<9F(-<p)/<9<p -F(-<p)^ p = V($)F(-f) 

• tachyon [26], p = V((j))/^\ - 2 p = -V{$)y/l - 2 

. quintom[27], p = (^ 2 -^ 2 )/2 + y(^^ 2 ) p = {$-$)/2-V{$ u <h) 



This schematic account of the possible types of matter/energy content in Cosmology 
does not make justice at all to the vast literature on the subject, but we must stop here 
and return to geometric aspects. 

Commonly, cosmological parameters are constrained by studying how distant light 
sources are seen from our detection devices, and this depends on the cosmological 
model. In an expanding FRW model light gets redshifted during its trip from the source 
to the observer, and the redshift depends on the amount of expansion occurred mean- 
while: 

a{tobs)/ a {temit) = l+Z = ^ bsl Kmit- (12) 

As light travels along null geodesies, one can compute the redshift experienced by a 
photon as it travels a given radial distance 

dt 2 — a 2 dr 2 = implies dr = dz/H(z), (13) 

we have taken ciq = 1 and c = 1 in dr = cdz/H(z) and the metric is the FRW one. Those 
are basic ingredients for the construction of "distances" in cosmology (see [28] for a 
pedagogical acount of the topic). 

Now, there is a basic concept we need so as to make progress in the topic of distances 
in Cosmology: comoving coordinates. One of the niceties of General Relativity is it 
allows to formulate physical laws in whatever system of coordinates one may prefer. 
In Cosmology, those coordinates have become very popular because they make one's 
life a lot easier. Observers which see the Universe as isotropic have constant values of 
spatial coordinates when the comoving system is used. Non-comoving observers will see 
redshifts in some directions and blueshifts in others, and the time measured by comoving 
observers is cosmological time. 

One of the quantities required for our geometrical tests is the line of sight comoving 
distance. This is the comoving separation between two objects with the same angular 
location but different radial position (P\Pz in the picture). The comoving distance from 
us to an object at z 

D c = f Z dz/H(z), (14) 
Jo 

whereas from the former we can derive that the comoving distance between two objects 
separated Az is Dc = Az/H(z). The comoving distance is the proper distance between 
those objects divided by the ratio of the scale factor of the Universe at the epochs of 
emission of reception. 

Another relevant definition is that of (transverse) comoving distance and angular 
diameter distance. The comoving distance between two objects located at the same radial 
position but separated by and angle AO is D M A6 (P3P4 in Fig. 2(b)) where D M = S(r) is 
the transverse comoving distance. Closely related to the former we have the widely used 
angular diameter distance Da - It is defined as the ratio of an object's physical transverse 
size to its angular size (in radians): 

D A = D M /\+z. (15) 

Finally, let me present the luminosity distance which is the key to the extraction of 
information from supernovae data. Given a standard candle its bolometric (i.e. integrated 



over all frequencies) emitting power right at the position of the source is called the 
luminosity L. The total bolometric power per unit area at the detector is called the flux 
F. The quantities L and F are used to define the luminosity distance Di = ^JL/AkF . 
At the moment of detection photons are passing through a sphere of proper surface area 
4n(a(t t, s )S(r)) 2 . The flux is affected by redshift in two ways: the energy decreases by a 
factor (1 +z) and the arrival dates are reduced by a factor (1 +z) so the flux is (1 +z) 2 
times smaller than in a static universe, so finally 

D L = (\+z) f Z dz/H(z). (16) 
Jo 

Actually, there is yet one more definition of distance which is used in one of the tests to 
be discussed below, but I prefer to postpone its mention for now. 

SUPERNOVA AS DARK ENERGY PROBES 

Theory has played the preeminent role in the development of cosmology till just a few 
years ago, with Einsteins' theory of General Relativity being the most compelling frame- 
work for the study of the evolution and fate of the Universe on large scale (we exclude the 
origin of the Universe from the list as the necessity to account for quantum effects makes 
General Relativity insufficient). However the preminence of theory seems doomed, as 
of recent the situation has changed greatly with the rise of advanced technological re- 
sources. Observations allow backing up predictions of General Relativity precisely and 
also make researchers formulate new questions. 

Daily experience connects us with the attractive side of gravity. It keeps us attached to 
the ground, it allows for the fun in games involving balls, and it holds artificial satellites 
revolving around the Earth. Yet gravity has a repulsive side which fits in Einstein's 
theory of gravity, but it was for decades believed to be just a theoretical possibility. No 
wonder, though, as it manifests on cosmological scales only. The existence of a exotic 
component in the cosmic budget which makes the Universe accelerate was inferred from 
the observation of distant supernovae [8, 9] 

In order to understand the evidence found in those experiments and why it was so im- 
portant, some background material is required. Supernovae are spectacularly luminous 
objects arising from explosions of massive supergiant stars. In the explosion a lot or all 
of the star's material is spelled out at a velocity of up to a tenth the speed of light. These 
phenomena have long intrigued astronomers, for instance Chinese records date back to 
AD 185, whereas first european record dates back to AD 1006 [30]. Supernovae a are 
classified according to the the shape of their light curves and the nature of their spectra 
[29], and they are fantastic for Cosmology as they may shine with the brightness of one 
thousand million suns and release a total of 10 44 joules (you could use it to provide the 
average USA consumption for 10 10 trillion years). Unfortunately supernovae are rare, in 
a given galaxy they only occur twice every thousand years, so there not as many of them 
available as cosmologists would wish. 

Let us now present some basic facts about supernovae. Type la supernovae (SNe 
la) are typically 6 times brighter than other supernovae, so predictably those are the 
most frequent supernovae in high redshift surveys. These explosions occur in binary star 



systems formed by a carbon-oxygen white dwarf and a companion star. The white dwarf 
accretes mass from the companion and when the Chandrasekhar limit is reached (10 14 
solar masses), the nucleus gets fused suddenly and the star explodes getting completely 
disrupted. 

An important part of the game are the different tools to determine the distance to a 
given supernovae (or to an astronomical object in general). The magnitude of a star is 
defined through its flux as m = — 2.5 log 10 F + const, so the brighter the star the lower 
the magnitude (there are plenty of places where this is explained, one is [32]. The 
photometric zero point [31] is related to Vega so the magnitude difference between two 
stars is 

mi -m 2 = -2.51og 10 Fi/F 2 . (17) 

The mnemotecnic rule is that a flux ratio of 100 gives a magnitude difference of 5. As 
flux F is related to luminosity distance through F = L/AtzD\, for two objects of the 
same intrinsic luminosity mi—m 2 = 51og 10 D L1 /D L2 - The absolute magnitude of a star 
is denoted by M and is defined as its apparent magnitude at a distance of 10 parsecs 
(32.6 light years) so m — M = 5\og l0 (Di/l0pc) and the quantity m — M is called the 
distance modulus. 

Another important aspect of the problem is to what extent supernovae can be consid- 
ered as standard candles. SNe la are not truly standard candles as they do not display 
the same luminosity at maximum and their appearance is not uniform. However, on the 
safe side we can say supernovae are standarizable as they can be brought in line with 
each other by some corrections: these are the "stretch factor correction" [33] and the 
'K-correction" [34], which are respectively a stretching or a contraction of the timescale 
of the event and a correction to compensate for the slight differences in the part of the 
spectrum observed by filters used to observe high and low redshift supernovae. 

Now let us enter the core of this section which is how to set up constraints on the 
parameters of the Universe using supernovae. As we have just mentioned, supernovae 
are not standard candles, but their luminosity curves can be unified to obtain a template 
value of the absolute magnitude. Photometry gives us the apparent magnitude m and our 
massaging of the light curves allows fixing M. This is only half of the story, though, 
because we also need to associate a recession velocity (or redshift z to each supernova), 
get rid of nuisance parameters (if possible), and formulate one's preferred model and test 
it. These were basically the steps followed in the pioneering works of 1998, and these 
are (more or less) the same steps everyone else in this factory keeps doing. 

In what recession velocities are concerned it must be kept in mind that astronomical 
objects have characteristic spectral features due to the emission of specific wavelengths 
of light from atoms or molecules. Templates exist which can be used to determine which 
kind of object one is observing one basically compares different spectra till a best match 
is found and which is the object's redshift. Supernovae are assigned the redshift of their 
host galaxy; Fig. 3 (a) is a courtesy of SDSS l . 
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(LAMOST), Los Alamos National Laboratory, the Max-Planck-Inslilute for Astronomy (MPIA), the Max-Planck-Inslilute for Astrophysics (MPA). New Mexico Slate University, Ohio Slate University, University 
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FIGURE 3. (a) A galaxy spectrum at four different redshifts (0.0, 0.05, 0.10, 0.15, 0.20). (b) Mock 
Hubble diagram 



The trick now is to draw the Hubble diagram of your supernovae sample(s) (see 
Fig. 3(b)). Recall that the larger m the fainter the object. At low redshifts (z < 0.1) 
one has a linear law m — M = \x ~ 5log\oz + jUo- But at high redshifts it is no longer 
valid and m — M = /! ~ 5logiQdi{z) + Mo an d the role of the matter and dark energy 

content become substantial. As one can deduce d^z) = z + (1 — qo)z 2 /2 + , and 

observations indicate cIl(z) > z at high redshifts, this is evidence indicating go < 
i.e. a currently accelerating universe. As qo depends on the cosmological model (e.g. 
qo = £l m j1 — Q.\ in LCDM) it gives an observational way to constrain the universe. 
On the other hand, by definition di(z) = (1 +z) J z Ho/H(z), from where it follows that 
what the observations are telling us is that H(z) > Hq in the past, which is another way 
of saying the rate of expansion has increased. 

Those best fits are obtained by minimizing the quantity 



the cfy j being the measurement variances [36]. The nuisance (statistically non- 
important) parameter /lo encodes the Hubble parameter and the absolute magnitude 
M and has to be marginalized over [35], so one will actually use be working with the 
quantity 



These % and o quantities deserve an explanation but I will postpone it for a few sections 
yet. 

Finally, in case you want to undertake studies of this sort by yourself you may like 
to know that at the time of writing this contribution the Davis et al. 2007 dataset [36] is 
one of the latest supernovae catalogs to be completed. It consists of 192 SNe classified 
as type la up to a redshift of z = 1.755. It is formed by 60 ESSENCE supernovae [37], 
57 Supernova Legacy Survey supernovae [38], 45 nearby supernovae [39, 40, 41] and 
30 Hubble Space Telescope (HST) supernovae [42]; and it is available at [43] Now, even 
Supernovae luminosities give the strongest evidence of the current acceleration of the 
universe, (and letting alone the problem with the physics of supernovae themselves), 
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(18) 





' ' ' I ' ■ ' ' I ' ' ' ' I ' ' ■ ' I ' ' ' 1 I ■ ' ' ' ' 
12 3 



FIGURE 4. 68% and 95% confidence regions of Q. m and Ha from various current measurements, and 
expected confidence region from the SNAP supernova program. Image courtesy of SNAP. 



there is the problem of a certain degeneracy in the test (see Fig. 4) so supernovae do 
not provide a good individual estimate of the cosmological parameters. We will address 
some of those tests in the next sections. 

HUBBLE PARAMETER FROM STELLAR AGES 

There is a feverish activity in constraining parameters of the Universe. We are focusing 
in these lectures in the tests based solely on geometry. The first goal of these researches 
is determining the current value of the equation of state parameter w. The next goal is to 
determine its evolution, i.e. to draw its redshift history. The supernovae luminosity test 
is somewhat compromised by its integral nature: cIl{z) = Hq{\ +z) /q(1 + z')(dt/dz'). 
As an alternative, measures of the integrand of the latter have been proposed in [44] as a 
way to improve sensitivity to w(z). That test relies on the availability of a cosmic clock 
to measure how the age of the Universe varies with redshift. Such a clock is provided by 
spectroscopy dating of galaxy ages. The idea is to obtain the age difference At between 
two-pas sively evolving galaxies born approximately at the same time but with a small 
redshift separation Az. One will then use At/Az to infer dt/dz, as it is directly related to 
the Hubble parameter H(z) = — (l+z)~ dz/dt, which is just the inverse of the integrand 
in the luminosity distance formula. 

However, in Astrophysics important open questions about galaxy formation and evo- 
lution remain. How did a homogeneous universe become a clumpy one? Inflation drives 
a major early homogeneization of the Universe but it also provides the primordial fluc- 
tuations which give rise to structure. How were galaxies born? How do they change over 



time? Progress in understanding these questions will help making the most out of this 
test. 

Hierarchically superclusters, clusters, galaxies, star clusters and stars can be distin- 
guished in what we are concerned in this section. Clusters may contain thousands of 
galaxies, and are often associated in larger groups called superclusters. Globular clus- 
ters are spherical collections of stars orbiting around a galactic core (so they are its 
satellites) are very tightly bound by gravity, and the oldest of them contain the oldest 
stellar populations known. They have been estimated to have formed at z > 6. It is these 
objects one resorts to in order to determine H(z). 

TABLE 1. H(z) data (in units of kms^Mpc - 1 ) 

z 0.09 0.17 0.27 0.4 0.88 1.3 1.43 1.53 1.75 

H(z) 69 83 70 87 117 168 177 140 202 
a 12 8.3 14 17.4 23.4 13.4 14.2 14 40.4 



On the other hand, most galaxies in cluster are elliptical galaxies, and of course this 
stands for the oldest galaxies in clusters too. This allows for some simplifications in their 
spectroscopic analysis. The derivative of cosmic time with respect to redshift is inferred 
from the aging of stellar populations in galaxies. Even though birth rate of galaxies 
is very high at high redshifts, examples of passively evolving (typically old and red) 
galaxies are known. 

In [45] three such samples were used: the first one was a sample of old red galaxies 
from the Gemini Deep Deep Survey, the second one was the so called Treu sample 
(made of field early-type) galaxies, and finally they used two radio galaxies (53W091 
and 53W069). A collection of 32 galaxies was completed this way and their ages were 
estimated using SPEED stellar population models [46] which were confronted with 
estimates by the GDDS collaboration [47] and were shown to give a good agreement. 

The next task was deriving differential ages, and this was done in various steps. 
Firstly, all galaxies within Az = 0.03 of each other were grouped together. This allowed 
estimating the age of the Universe at a given redshift. The redshift interval is small 
enough to exclude galaxies which have undergone a significant age evolution, but large 
enough as for the bins to be made of more than one galaxy. Secondly, age differences 
were calculated by comparing bins with 0. 1 < Az < 0. 15. The lower bound results in an 
age evolution larger that the error in age determination; Thus, one achieves a robust age 
determination. Finally, the value of H(z) was derived using the convenient expression 
given above, and then Table 1 was constructed. 



COSMIC MICROWAVE BACKGROUND 

The first announcement of the discovery evidence of the cosmic microwave background 
(CMB) dates back 42 years. The radio astronomers Penzias and Wilson did the the 
discovery but did not realize what it was due to. In broad terms they detected an 
extraterrestial isotropic excess noise using an Bell Labs antenna sited at Holmdel (New 
Jersey) which after years of being used for telecommunications had by 1962 been freed 
up for pure research Dicke's group at Princeton (formed by Peebles, Roll, Wilkinson 
and Dicke himself) gave the theoretical interpretation of the Penzias and Wilson result 



as a prediction of the big bang models. Two papers [4, 5] were sent jointly and appeared 
in the 142th volume of the Astrophysical Journal, the Penzias and Wilson paper was 
humbler, but they got rewarded years later (in 1978) with the Noble Prize. Theoretical 
advances and predictions on the topic had been done earlier by Gamow, Alpher and 
Herman, and although the Princeton group had drawn their conclusions independently, 
one of his members, Peebles, recognized the contribution of those authors in his classical 
textbook [48]. Just a few months after the publication of the mentioned two papers, 
Roll and Wilkinson [6] confirmed in another work the thermal nature of the spectrum 
of the radiation (in concordance with the big bang model prediction) and estimated its 
temperature to be about 3.0 ± 0.5 AT. 

The next major breakthrough in this topic was the discovery of the CMB anisotropics 
by the COBE satellite, and as a consequence Cosmology got two more Nobel laureates in 
2006: Smooth and Mather. According to the calculations carried out by Harrison [49], 
Peebles and Yu [50], and independently by Zeldovich [51], CMB anisotropics are the 
consequence of the predictable primeval inhomegeneities of amplitude 10~ 5 . 

The CMB power spectrum, through its pattern of peaks (see Fig. 5(a)), tell us that the 
Universe is very close to spatially flat (at a high degree of accuracy) (see for instance 
[52] among the also vast amount of references mentioning this). The CMB also inform 
us that inflation was the major component of cosmic structure formation as opposed to 
cosmic strings (for a short history of the discovery of the first peak see [53]) 

You may already be convinced that CMB physics is capital. On the one hand no other 
cosmological probe has beaten it so far in the amount of relevant information it provides, 
and on the other hand it supports the big bang model and links Cosmology and particle 
physics as the observed abundances of light elements determined from measurements 
of the fluctuations of the cosmic background radiation temperature are in agreement 
with big bang nucleosynthesis [54], being based on the well-trusted Standard Model of 
Particle physics. But the list of important pieces of information it gives us does not end 
there as it provides traces of very weak primordial fluctuations which may have seeded 
the large scale structure observed today, it informs us of the geometry of the universe 
(first peak), it provides evidence for baryonic dark matter (second peak), it constrains 
the amount of dark matter (third peak) [55] and finally, a fact which is very important to 
us is that constraints on dark energy models can be refined using CMB data. 

At this stage is it convenient to make a thermal history sketch (as illustrated in Fig. 
5(b)). After the big bang the Universe cools down due to adiabatic expansion and its goes 
through these stages: the energies of particles decrease and several phase transitions 
occur, massive particles become non-relativistic as temperature decreases below their 
rest mass, the rates of production of particle-antiparticle decrease and anhilitation leaves 
asymmetric populations, nuclei of light elements form (z — 10 4 ) and the Universe 
becomes a soup of nucleons, photons and free electrons (the CMB forms at z — 10 7 ). 

As the density of relativistic particles decreases faster (proportional to a~ 4 ) than that 
of non-relativistic ones (proportional to a~ 3 ) at some stage one gets matter-radiation 
equality. Later on, thermal ionization ceases to be effective, atoms form (z — 10 4 ) 
electrons get bound to nucleons (recombination) and photons decouple and travel freely 
to us since [56, 57] 

Just to give a few more details about the formation of the CMB background let us 
mention that it formed during the epoch 10 10 > z > 10 7 due to photon creation processes 




(bremsstrahlung and double Compton scattering) [58]. At later stages the dominant 
coupling effect turned to be Compton/Thompson scattering but it did not destroy the 
thermal nature of the spectrum. 

The basic observable in CMB physics is temperature fluctuation AT /T (a field de- 
pending on direction on the sky). Initial density fluctuations were tiny, with an ampli- 
tude of about 10 so their evolution is linear, i.e. Fourier modes evolve independently. 
These primordial irregularities give rise to potential wells and hills according to which 
the baryons and the photons get accommodated. Photons, climb up (down) potential 
wells (hills) in an unimpeded trip to us after getting released at the decoupling redshift 
Zdec ~ 1100 and give us a snap-shot of the early universe in the form of cold and red 
spots. This is a primary anisotropy that gets imprinted on the CMB at last scattering 
responsible for the characteristic large scale anisotropy called Sachs-Wolfe plateau [59] 
(this anisotropy is associated with fluctuations with period so large they have had no 
time to perform a complete oscillation by the recombination time [60]). 

Let us now discuss were the peaks and troughs come from. Photon pressure supplies 
force: Vp r = Vp y /3, py ~ T 4 , and the spatial variations of density themselves are 
also a source of changes in the gravitational potential. Potential and pressure gradients 
compete between them and induce acoustic oscillations in the fluid which result in 

temperature oscillations AT ~ §py^ ~ A(k)cos(kc s t), (harmonic wave), with c s the 
speed of sound. This photons do not find opposition either between their departure and 
their arrival to their destination (our observing devices). Photons released at maximum 
compression give blue spots as they have to climb up the potential, so they appear blue 
to us. In contrast, those caught at maximum rarefaction form red spots. In summary, the 
pattern of spots associated with long wavelength perturbations gets enriched by those 
coming from perturbations with smaller wavelengths. 

That of course, is not the end of the explanation, the angular power spectrum is the 
function describing how the variance of the fluctuations depends on the angular scale (its 
definition rests on the assumption of Gaussianity and randomness of the fluctuations). 
Fluctuation extrema correspond to peaks in power (even ones correspond to compression 
and odd ones correspond to rarefaction). Randomness of the fluctuations implies we 
have them in all wavelenghts, and clearly and when we compare a fluctuation of a given 



wavelength with another one with half of that wavelength it is obvious that the time the 
second takes to complete an oscillation is half of the time the first one takes. Thus, if by 
recombination there is a mode which has completed half an oscillation, there will modes 
which will have completed 2, 3, 4, . . . half oscillations. 

You may have noticed, as well, that the spectrum presents a modulation, it is due to 
different effects (see [61] for an authorized discussion). Baryons being massive have 
preference for the throughs and the compression gets enhanced, thus the amplitude of 
odds peaks gets larger. Another effect of baryons is making the oscillations slower so 
the peaks are pushed to higher multipoles. On the other hand, photons are also important 
in the modulation as in the radiation dominated phase. During radiation domination 
most gravity comes from photons, they are the major source of gravitational potential 
so as pressure redistributes the photons, gravitational potentials get washed away. This 
effect is absent when the Universe becomes matter dominated, so it is exclusive to high 
frequency modes, whereas low frequency modes do not start to oscillate till matter 
domination has begun. Finally, there is also a strong damping effect due to slight 
imperfections in the fluid associated with shear viscosity and heat conduction (Silk 
damping [62]). 

Interestingly, the CMB provides a dark energy test which requires considering only 
the background geometry, in contrast to other uses which rest on the powerful but chal- 
lenging study of perturbations. This is the CMB shift test, which estimates how a physi- 
cal length in the primordial universe appears to us. By comparing the (non-Euclidianity) 
effects on geometry of different matter/contents inference about the likeliest ones can be 
done. The basic quantity to consider is that in an arbitrary model the position of the first 
peak in the temperature spectrum is l\ = kDa{z kc ) / r s {z K c), where the quantity r s {z KC ) 
represents the last scattering sound horizon scale i.e. the distance sound may travel be- 
fore the recombination epoch. 

The CMB shift gives the first CMB peak position ratio between a model one wants to 
test (unprimed model) and a reference Einstein-de Sitter model (primed model) 

R = 2h/l x 

[63] It is considered as a robust test as it does not depend on the parameter Hq (the 
Hubble factor today), which is not constrained by the other tests. Considering the speed 
of sound c s is constant and using the approximations 

D' A (z re c ) ~ 2ca rec /Hq r s (z rec ) Rj 2c s a r lc / 'Hq \/^m , (19) 

one finally arrives at 

R^Hq^Q^ dz/H(z). (20) 
Jo 

You can go by the value i?(1089) = 1.71 ±0.03 calculated from WMAP3 data. 

In a fashion similar to the other tests one will have to construct a x 2 function using 
the latter observational value and the theoretical function for the cosmological model to 
be constrained. 

In the next section we are going to consider a test which also stems from early universe 
physics. 



FIGURE 6. Two stages in the evolution of a primordial density perturbation. Baryon density is on the 
left panels, photon density middle ones and mass profile in the right panel. The upper image is for a 
situation not much later than the decoupling time, the last one correspond to much later epoch. Images 
courtesy of Martin White. 

BARYON ACOUSTIC OSCILLATIONS 

Diagnosing the mystifying new physics causing acceleration requires more precise mea- 
surements of cosmological distance scales. In recent years a new geometrical constraint 
of dark energy has emerged which relies on traces left by early universe sound waves in 
the galaxy distribution. This test is somehow in its infancy but it is very promising. As 
it often in science, quite a few years have elapsed form prediction of the effect [64, 65] 
till its detection [66]. 

Let me move on now to a physical description of the phenomenon based on the one 
found in [67] (see as well [68, 69]). The early universe was composed of a plasma of 
energetic photons and ionized hydrogen (protons and electron) in addition to other trace 
elements. Imagine now a single perturbation in the form of a excess of matter. Pressure is 
very high and ejects the baryon-photon outward at relativistic velocities. At first photon 
and baryons go by the hand (the speed of the radius of the shell being larger than half the 
speed of light). At recombination photons decouple and stream away whereas the baryon 
peak gets frozen; the reason for this is they are cold dark matter and therefore they have 
no intrinsic motion (no pressure) unlike the previous stages where they were coupled 
to the photons (in this early situation baryons were subject to radiation pressure). The 
photon distributions become more and more homogeneous, while the baryon overdensity 
does not disappear. Finally, the initial large gravitational potential well begins to attract 
matter back into it and a second overdense region appears at the center. Two stages in 
the evolution of a single perturbation are illustrated in Fig. 6 

The most important fact of the phenomenon is it leaves traceable effects on large-scale 
structure. Initial perturbations would have produced wavy excitations emanating from all 
points, and the subsequent evolution would resemble what happens when dropping rocks 
in a pond. The plasma shells would have today 150 Mpc (or 500 million light years) 
and galaxy formation in the locus of those shells would be likelier. In fact, there is a 
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FIGURE 7. Correlation function versus redshift. Image borrowed from [66]. 

correlation between galaxies in the shells and their centers which result in the detectable 
effect that galaxies are more likely to be separated by that distance than by larger or 
smaller ones. 

In consequence, the large scale two-point correlation function presents a clear peak 
at 100/? 1 Mpc. The two point correlation B, function (see Fig. 7) controls the joint 
probability of finding two galaxies centered within the volume elements dV\ and dVj at 
separation r (the galaxy number density is n and ro is the characteristic clustering length): 
dP = n 2 [l + %(r/ro)]dVidV2 The derivation of the scale of the shell r pea ^ involves 
postulating the form of H(z), and as the value of r pea ^ predicted by CMB physics is 
known one can vary the parameters of the model to find the best fit. The conversion 
between distances and redshifts can be done using the dilation scale D v (z). This quantity 
is an "isotropic" distance definition encoding the effects of expansion along the line of 
sight and along the transversal direction to it (see Fig. 2(b)): 

D v (z) = [((P 3 P4)/A0) 2 (PiP 2 )] 1/3 = (D M (z) 2 ADc) l/3 (21) 

The newest results are those given in [70]. Three galaxy catalogs were used: a main 
galaxy catalog combining 2dFGRS and SDSS main galaxies (this catalog fixes Dy (0.2)), 
a catalog of SDSS large red galaxies (this catalog fixes Dy (0.35)), and a combi- 
nation of the two. This combined sample gives r s /Dy (0.35) = 0.1094 ± 0.0033 and 
r s /D v (0.35) = 0. 1980 ± 0.0058. This data must be inserted in the adequate x 2 function 
and it will only remain to choose what model to test. 

COSMOLOGIES WITH LATE-TIME ACCELERATION 

There are, as you can imagine, a lot of cosmological models in the market with more or 
less physical foundation you could try and apply the dark energy tests presented here. 
Basically, the lesson to be extracted from the previous discussion is that in principle just 



needs to postulate the functional form of H(z). In what follows I wish to make an a 
outline of a few ot these models. This is, of course, a biased list of parameterizations as 
it is based on my own preferences, and on the other hand there are authorized reviews 
you can head to for a more exhaustive revision. 

In general, a given general relativistic model of dark energy with equation of state 
Wde (z) leads to 

// 2 (z)/// 2 = Q m (l+z) 3 + (l-Q m )exp 3 [\l+w de (x))/(l+x)dx . (22) 

To get the latter a universe containing dark energy and dust (dark matter and baryons) 
has been assumed. 

The first model I want to consider here is the LCDM model. It follows from the 
choice w c ie(z) = — 1, i.e. dark energy is a cosmological constant. Thus 



h 2 { z )/hI = a m {\ +z) 3 + i -a„ 



(23) 



with Q. m a free parameter. LCDM is consistent with all data, but there is a theoretical 
problem in explaining the observed value of the cosmological constant A. 

The second case of this list is the DGP model It was proposed by Deffayet [71] 
the inspiration being [72] It represents a simple alternative to the standard LCDM 
Cosmology, with the same number of parameters. Late-time acceleration in this model 
is due to an infrared modification of gravity (no dark energy needed). Explicitly one has 



h(z)/h = (i - a m )/2 + ^(i-a m ) 2 /4+a m (i+z) 3 . 



(24) 



Let me now tell you about the QCDM model. It is the simplest generalization of the 
LCDM model. It consists in taking a constant value of Wd e = w different in general from 
— 1, so 



H 2 (z)/Hi = Q ra (l + z ) 3 + (1 - a m )(l +zf l+ 



w) 



(25) 



Even though it is a rather simple model, it can be useful for detecting at a first instance 
preference for evolutionary dark energy. 

Another two parameter model of interest is the LDGP model [73]. The DGP model 
has actually two separate branches They reflect the two ways to embed the 4D brane 
universe in the 5D bulk spacetime. We have loosely called before DGP before to the self- 
accelerating branch (this is common practice in the literature) but I insist there is another 
branch which is physically interesting too. The LDGP model precisely represents this 
branch, i.e. the non self-accelerating one, with a cosmological constant included which 
is required for acceleration. Interestingly, higher values of the cosmological constant are 
allowed as compared to LCDM (due to an screening effect). In a explicit way 



H(z)/H = ^Cl m (l+z) 3 + l-a m + 2v^T + 



(26) 



and the parameter Q. rc encodes the so called crossover scale signaling the transition from 
the general relativistic to the modified gravity regime 



The last but one case I wish to address is Chevallier-Polarski-Linder Ansatz [74]. 
This is a widespread generalization of QCDM for which 

w de {z) =WQ + W\(\ - 1/(1 +z)). 

Therefore 

H 2 (z)/Hi = a m (l +z) 3 + (1 -n m )(l +z ) 3 ( 1 +-o+-i) e -3- lZ /i+z. (27) 

Two nice properties of it stand out, firstly the model displays finiteness of H(z) at high 
redshifts, secondly it admits a simple physical interpretation as w\ is a measure of the 
scalar field potential slow roll factor V'/V in a quintessence picture [76] 

Finally, to put and end on this list, I would like to consider the QDGP model 
[75]. This is a generalization of LDGP: the cosmological constant is replaced by dark 
energy with a constant equation of state w. The modified 4D Friedman equation has the 
following form: 

H(z)/H = Jajl+z) 3 + (1 -£2 m + 2v^)(l +z) 3 ^ +£l rc - v^-(28) 

A remark about three of the models considered is in order. Models with late-time 
acceleration due to infrared modifications of gravity are called dark gravity models, and 
DGP, LDGP and QDGP are of that sort. An effective dark energy equation of state w e ff 
can be deduced by imposing p e ff + 3H(1 + w e ff)p e ff = 0. Calculating this in terms of 
redshift for the three models considered here can be a nice exercise for the readers. 

I finish here this section as a turn of subject is required, as I feel this contribution 
would not be sort of self-contained if the issue of the statistical treatment of the tests 
considered was not covered. 

STATISTICAL INFERERENCE 

Science, and ergo Astronomy, is about making decisions. All scientific activities (de- 
sign of experiments, design and building of instruments, data collection, reduction and 
interpretation, ...) rely on decisions. In turn, decisions are made by comparison, and 
comparison requires a statistics (a summarized description) of the available data. The 
term "statistics" means rigorously a quantity giving broad representation of the data, but 
the same term is also loosely used when referring to "statistical inference". Moreover, 
as in science making decisions involves measurements and derivation of values of pa- 
rameters, it is necessary to assess the degree of belief of the true value of the parameter 
being measured/derived. 

Astronomy/Cosmology being about the study of the Universe makes them singular in 
what concerns decision making. The Universe is not an experiment one can rerun, and to 
make things worse the typical objects studied are very distant. In consequence, typically 
one will have very poor knowledge of the distribution behind the variables subject to 
measurement. Yet the difficulties must be overcome as statistics enters all five stages in 
the loop characterizing every experiment: observe-reduce-analyse-infer-cogitate. 



I am going to introduce now some basic notions, definitions and notation to be used 
here after. "The readers interested in learning more about most of the topics addressed 
in this section can head to [77], but I also recommend to you other sources such as 
[78,79,81]. 

Probability will be a numerical account of our strength of belief. Let us enunciate now 
Kolmogorov axioms: for any random event A one has < prob(A) < 1, if and event A is 
certain then prob(A) = 1, if event A and B are exclusive, then prob(AorB) = prob(A) + 
prob(B). These axioms are basically all is needed to develop entirely the mathematical 
probability theory (by this one means the recipes to manipulate probabilities once they 
have been specified). Probabilities are not inherent properties of physical problems, that 
means they do not make sense on their own, they are just a reflection of the knowledge 
we have. 

If the knowledge about even A does not affect the probability of event B, those 
events are independent, i.e. prob(AandB) = prob(A)prob(B). If A and B are 
not independent events then it is convenient to know conditional probabilities: 
prob{A\B) = prob{AandB)/prob{B) and prob{B\A) = prob(AandB)/prob(A), 
which can be derived from Kolmogorov axioms. If A and B are independent, then 
prob(A\B) = prob(A) prob(B\A) = prob(B) In addition, if event B comes is dif- 
ferent flavours (5i,52,#3, •••)> then p(A) = Y*iProb(A\Bi)prob(Bi). Moreover, if A 
is a parameter of insterest but the B t are not, knowledge of prob(Bi) allows getting 
rid of those nuisance parameters by summation/integration, and the process is called 
marginalization. 

Having presented this basics let me guide you into the realm of Bayesian infer- 
ence By equating prob(AandB) and pwb(BandA) one gets the identity prob(A\B) = 
prob(B\A)prob(B)/ prob(A) known as Bayes theorem. 

There is a state of belief, the prior prob(B), before the data, the event A, are col- 
lected, and experience modifies this knowledge. Experience is encoded in the likelihood, 
prob(A\B). The state of belief at the end of the process (analysis of the data) is repre- 
sented by the posterior prob(B\A). Mathematically naive though it is, this theorem is, a 
regards interpretation, very powerful, but not devoid of controversy. I provide here only 
a short account of Bayesian statistics, if you feel like learning more about it, head to this 
site [80]. 

For a better understanding of the different between the Bayesian and frequentist 
approaches, consider now taking out balls from a box (you cannot see the inside and 
only know they are either white or red). Imagine you extract several times three balls and 
put them back into the box. Would you care of the probability of taking out two whites 
and one red? Would you rather not be more interested in saying something about the 
contents of the box? The second question seems the more natural one. This is the kind of 
question Bayes theorem allows to answer. Bayesians care for probabilities of hypotheses 
given data, whereas frequentists are concerned about the probability of hypothetical data 
assuming the truth of some hypothesis. 

A useful concept in connection with all this is that of probability density functions. 
If x is a random real variable expressing a result, the probability of getting a number 
near x is p(x) = pwb(x)8x, and p(x) is the probability density function distribution 
of x (probability density functions are also called probability distributions). Probability 



density functions satisfy these properties: prob(a <x<b)= f% p(x)dx, J^ oa p(x)dx = 1, 
p(x) is single- valued and non-negative for all real x. Probability density functions are 
usually quantified by the position of the center or mean /i = J co 00 xp(x)dx, the spread 

o 2 = p oa (x-n) 2 p(x)dx. 

Among them one stands out notoriously; it is the Gaussian probability density func- 
tion, which is an ubiquitous distribution in Physics, thanks to the central limit theorem, 
which makes it work is most situations: broadly speaking a little averaging makes any 
distribution converge to the Gaussian one, explicitly, p(x) = e~( x ~^ l G / \J2%<5. 

We are going to make use of it for parameter estimation and model selection, which 
is at the core of many investigations in Cosmology these days. Consider we have N 
observational data for some physical quantity f th of interest. We set {f° bs } = {dj}. 

Assume f h depends on parameters {Q t } and consider it in the context of a given model 
In addition, regard the errors as Gaussian distributed so the likelihood reads 

&{{d J }\{6 l },jr) « e -* 2 « e <» with X 2 ({ei}) = £ [fj bs -f h ({ei})} 2 /aj. (29) 

7=1 

For uncorrected observational datasets 

^({jj 1 V • • • n {4 m) )l W^) ~ 

The probability density function p({6j}\{dj},^) of the parameters to have values 
{6i} under the assumption that the true model is M and provided that the available 
observational data are {dj} is given by Bayes theorem. In terms of probability density 
functions 

Here p({6i}\{dj},^) is called the posterior probability density function, and best fit 
values of the parameters are estimated by maximizing it, whereas 7t({0i},^) is called 
the prior probability density function. Choice of prior is subjective but compulsory in 
the Bayesian approach; remember the prior pdf encodes all previous knowledge about 
the parameters before the observational data have been collected. 

The first step to estimate parameters in the Bayesian framework is maximizing the 
posterior p({9i}\{dj},\J%). The second step is construction credible intervals It is 
convenient to simplify our notation, so p{{&i}\{dj} , = p(0i, ■ ■ ■ , 6 n ). The marginal 
probability density function on Q t is p{0i) = J p(0i, . . . , 9 n )dd\ . . .dOi \ dO i+ i . . .dd n . 
Of course, if the model is genuinely uniparametric no marginalization is required. It is 
also convenient to define parameters 9u and 6i u satisfying p(6u) ~ and p(0j u ) ~ The 
credible intervals, under the hypothesis that the marginal pdf is approximately Gaussian, 
are constructed from the median and errors. The 68% percent credible intervals on the 
parameter 0, are given as ; = x~t* 

The median x is calculated from f£ u p(0i)d6i = 0.5 x / p(6j)d0i The lower error y 
is calculated from fg~ y p(6i)d0i = ((1 — 0.68)/2) x J p{9i)d9i The upper error z is 
calculated from f°™ z p(0i)d0i = ((1 -0.68)/2) x / p(0;)</$ 



Credible contours are a popular/illustrating construction and they are worth a mention. 
By marginalization with respect to all parameters but two one gets 



P (e i ,e i+1 )=p(e i ) = J p(Oi,...,o n )do l ...do i - l do i+2 ...do n 

Again, obviously, marginalization is not necessary if the model is genuinely biparamet- 
ric. An effective x 2 can be defined in the form — 21og(p(0,, 0;+i)) = Xeff(®ii ^+1) /2 
(up to a constant), so the best fit responds to the maximum posterior criterion. The prob- 
ability that for some parameters other than those for the best fit (maximum likelihood) 
X 2 increases with respect to the best fit by an amount Ax 2 is 1 — T( 1 , \og(p b f/p))/T(l). 
By fixing the desired probability content, the latter becomes the implicit equation of the 
credible contours on the parameter space (popular choices in the literature are 68.3%, 
95.4% and 99.7%). 

Now let me comment about model selection. Bayesians use a estimator to select 
models which informs about how well the parameters of the model fit the data. It does 
not rely exclusively on the best-fitting parameters of the model, but rather it involves 
an averaging over all the parameter values that were theoretically plausible before the 
measurement ever took place [81]. Bayes evidence is 

g(JK) = p({dj}\^) = J Tt({di},JZ)^{{dj}\{d i \,Jt)dd l ...dd n , (32) 

so it is the probability of data {dj} given the model M Here 7l({6i},^) is the model's 
prior on the set of parameters, normalized to unity (i.e. / n({6i}, ^)d0\ . . .d6 n = 1.) 
Using the popular top-hat prior n(6i) = (0 jm ax — fy'min) -1 Bayes evidence is rewritten as 

(/"#lmax /"(Wax \ / /"^lmax f&nmix \ 

/ .../ j?(o h ...,o n )do l ...o n )/l ... dd x ...de n \. 
•^lmin ^fyrniin / V 1 ^ ^lmin •' fyimin / 

Finally, preference of model ^ over model <Mj given {<4} is estimated by 

p{JZi\{d k })/p{JZj\{d k }) = Si^Kii^/Sji^^Kji^j) 

The Bayes factor 5,j for any two models and is 

B ij = g i {Jt i )/gj(JZ j ) 

so if as in usual practice one assumes no prior preference of one model over the other, 
that is, Assuming = nj(^j) = 1/2 (no a priori preference) then 

p{.£i\{d k })/p{.£j\{d k }) =Bjj. 

The most popular key to interpreting Bayes factors is Jeffreys scale [82]: if In(Bij) < 1, 
then the evidence against is not significant, if 1 < In(Bij) < 2.5, then the evidence 
against is substantial, if 2.5 < In(Bij) < 5, then the evidence in favor of ^ is strong, 
if 5 < In(Bij), then the evidence in favor of M\ is decisive. 

There are of course many more things which could be mentioned about this topic, but 
I hope this condensed introduction will help you start crossing swords with the art of 
constraining dark energy models using geometrical tests like the ones discussed here or 
upcoming ones. 



CONCLUSIONS 



In this contribution I have tried to review some of the background on geometrical tests 
of dark energy models. A historical review of the development of Cosmology has been 
followed by an account of the importance of the discovery of the late-time acceleration, 
which is commonly attributed to the existence of an exotic component in the cosmic 
budget. Other possible explanations have been attempted, which required being open 
minded enough to admit the existence of extra dimensions. 

Examples of both conventional (if that adjective can be used) and extradimensional 
models have been mentioned here with respect to their geometrical features. All these 
models can be subject to different observational tests which only require postulating 
a parametrization of the Hubble factor of the Universe or quantities derived from it. 
Here I have been concerned with four tests only: luminosity of supernovae, direct 
H(z) measurements, the CMB shift and baryon acoustic oscillations. I cannot deny this 
selection is biased in the sense these are tests I have made research on, but the first and 
the last one are definitely very important and much effort is doing by research groups all 
over the world in finding the data their application requires. 

This contribution has also tried to maintain a pedagogical tone as it has been prepared 
for the Advanced Summer School 2007 organized by Cinvestav in Mexico DF. I just 
hope it will be useful to any reader which happens to come across it. 
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